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Q frared as a consequence of the amplification of superhorizon modes, the energy levels 
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I. INTRODUCTION 

It is expected that an observer, corresponding to an Unruh detector coupled to a scalar 
field, when freely falling in de Sitter spacejWill perceive radiation with a thermal spectrum of 
the de Sitter temperature Th = H/{2'k) 0,y], where H denotes the Hubble parameter. It is the 
purpose of this article, to clarify in what sense this result is universal to scalar fields of different 
couplings to the de Sitter background and how the detector apprehends the differences. 

Let us therefore refine what is meant by the observation of thermal radiation: At first order in 
perturbation theory, the detector response function is proportional to the Fourier transform of the 
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scalar propagator w.r.t. the proper time of the detector, and it describes how many particles are 
absorbed and emitted per unit time. When being in equilibrium with the de Sitter background, 
the energy levels of the detector are thermally populated according to the temperature T^- 

The easiest way to derive this is to note that the propagator for nonminimally coupled scalars 
has in the imaginary direction of its proper time t the periodicity t ^ t -\- 2-ni/TH 1^, As 
we will point out, this however does not completely characterize the response function of the 
detector, which describes the number of particles detected per unit time. The rate turns out to 
depend on the scalar mass and on its coupling to the curvature, as was first shown in Ref. 0], 
circumventing the use of the scalar propagator. 

The fact that in de Sitter space the invariance of the quantum vacuum becomes manifest when 
the scalar propagator only depends on the proper time seperation alon g a g eodesic has led to the 
practice of defining the de Sitter vacua through this quantity 0, 0, Ifli llfli However, in the 
case of a massless scalar, which is minimally coupled to the curvature, this leads to a problem 
since the propagator is infrared divergent 10|. We argue that, when regulated by a cutoff, this 
divergence gives rise to a contribution which is irrelevant to the total detector response. In 
addition, we compare with the response functions of a detector immersed in a thermal bath in 
Minkowski space and discuss the situation in de Sitter space-times with dimension other than 
four. 



II. UNRUH'S DETECTOR 

Unruh's detector P, Q, 0] corresponds to a heavy slowly moving particle along a trajectory 
X = x{t), which interacts with a scalar bath of particles as 

Cvnrnh = -hm{t)^{x{t)) , (1) 

where /i is a coupling constant, t is the proper time and rh{t) represents the quantum operator 
describing the monopole interaction of the detector with the scalar bath $ = Since the 

detector is assumed to be very heavy, it does not fluctuate in space, and hence the only time 
dependence is in rh{t). The response of Unruh's detector can be derived as follows. One assumes 
that the state of the detector is specified by a set of energy eigenstates, {l-E)}, and that each 
absorbtion of a scalar quantum elevates the energy of the detector by AE = E — Eq; the converse 
is true for each emission. 
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Consider now at first order perturbation theory the transition amphtude Ai from a state 

\Eo) (S> Iv^o) to a state {E\ (g) (v?! 

M{Eo-^ E;to,tf) = hmE,,E / ' rfti e-*(^-^°)*^(<^|$(a;(ti)) |(^o) , (2) 

J to 

where mEg^E = {E\rh{0)\Eo) (m(ti) = e*-'^*^m(0)e^*-'^*^) is defined such that it does not include 
stimulated emission. In order to take account of the influence of initial conditions, for simplicity 
we assume that the interaction m(ti) switches on at ti = as 0(ti — to); where 6(x) = 1 when 
a; > 0, and = when a; < 0. More generally, one would expect that the interaction switches 
on adiabatically. Since we are not interested in studying in detail how the response function 
depends on the initial conditions and on how the interaction turns on, this rough treatment 
should suffice. Upon squaring Ai and making use of orthonormality and completeness of the 
scalar eigenstates {Iv^)}, and taking to —oo, we get for the transition probability per unit 
proper time, 

. , W (^E^E-E„), (3) 

where 

dT{^E) ^ r dAte-'''^''HG<(x{t + At/2y,x{t-At/2)) (4) 

Clt J —oo ^ ' 

is the response function per unit proper time. Here zG^(xi;x2) = {ipQ\^{x2)^{xi)\ipQ) is the 
(positive frequency) coordinate space Wightman function, and we defined the time variables, 
t = (ti + ^2)72 and At = ti — t2. The expression (0)) is real by construction, which can be 
formally shown by taking the complex conjugate of the integrand in ^ and then making use of 
the hermiticity property of the Wightman functions, [iG'^'^{xi,X2)]* = iG^'^{x2, xi). We shall 
study the particle spectrum observed in an appropriate vacuum state \ipo) = |0), the precise 
nature of which will be specified later. It is generally accepted that it is most natural to consider 
(a freely-falling) Unruh detector moving on a geodesic, since in this case the detector does not 
see additional particles due to the Unruh effect, which would be present only because of the 
acceleration of the detector. 

A nice property of Unruh's detector is the separability of the transition probability © into a 
product of the selectivity function, which depends on the inner structure of the detector, which is 
completely specified by the operator m{t) and the coupling h, and the response function, which 
depends on the state of the scalar field only. It appears useful to rewrite Eq. (jH) in terms of the 



4 



Wigner function as 

f- = l0^s^G<{k^ = AE,kMt)), (5) 

where the Wigner function is defined as the Fourier transform w.r.t. the relative coordiante of 
the Wightman function, 

iG< {k, x) = j dW'-''iG< ix + r/2,x- r/2) , (6) 

such that r'^ = Xi — are proper times. This means that the response function of Unruh's 
detector is completely insensitive to particle momenta (which is consistent with the assumption 
that the detector must be very massive), and it measures (absorbs) scalar particles of all possible 
momenta k; likewise, it isotropically emits particles of all momenta. 



III. THERMAL RESPONSE 



In order to get an insight in how a certain distribution function of a scalar field is perceived by 
the detector, we consider the response function for the Bose-Einstein distribution. The Wightman 



function for a thermally excited scalar field has the well known form 
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zG<^{k) = 27rsign(A;°)5(A;2 + mj)-^^ , (7) 

where /? = 1/T denotes the inverse temperature, and the scalar mass. The response func- 
tion (0) for d dimensions in flat space-time is easily calculated, 

'^"■^f"" = ^^-gn(Aii;)e ({AEf - ml) {(AEf - mj) , (8) 

and explicitly, for d = 4, one finds 



d^,M^4(Ai^) _ ^ig^(^^)e ii^Er - ml) ^{AEy - ^ 

dt 2n e/5^^-1' ^' 

Hence, the response function of Unruh's detector to a scalar thermal state contains, apart from 
the Bose-Einstein distribution, an additional factor, which depends on the scalar particle mass 
and reduces to AE/{2tt) in the massless limit and d = 4. This is precisely the factor that arises 
for the conformal scalar vacuum (fT^ in de Sitter space P]. 



5 



IV. SCALAR FIELDS IN DE SITTER INFLATION 

The Lagrangean for a massive real scalar field in a curved space-time background is given by 



"^^^^ = -\^gg^''{d,^){d,^) - ^v^K + ^7^)$^ (10) 



where g = det[g^iy] denotes the determinant of the metric g^i,, g^" is the inverse of the metric, 
the scalar mass and TZ the curvature scalar. In conformal space-times, in which the metric 
is of the form 

gfiu = a^Vi^v , (11) 

where rj^y = diag(— 1, 1, 1, 1) is the fiat Minkowski metric and a = a{ri) the scale factor, the 
Lagrangean (|Tn|l reduces to 

= -laY''{d,<^){d.^) - la\ml + ^n)<l>'. (12) 

For example, in (a locally) de Sitter inflation a = -l/{Hr]) (r/ < 0), 7^ = 12H'^ (in D = 4 
space-time dimensions), while in radiation-matter era, a = a^r] + amv'^, where H denotes the 
Hubble parameter, rj conformal time, and are constants. 
The special cases of interest are: 

• conformally coupled massless scalar field (^ = 1/6, m,^ = 0), for which a conformally 
rescaled scalar 9? = a$ satisfies the simple differential equation (V denotes a spatial deriva- 
tive) 

{d',-V'Mx) = 0; (13) 

• and nearly minimally coupled light scalar (|^| ^ 1, rUfp <^ H), which obeys 

{d', - - ^0 + a\ml + e7^))^(x) = . (14) 
(in de Sitter inflation, —a~^<fa/d'rf = —2/iff = —2a?H'^)] 

• and the minimally coupled massless scalar (^ = 0, = 0), which satisfies the following 
differential equation. 
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A. Conformal vacuum in de Sitter space 

We now calculate the response function (0)) for a conformally coupled massless scalar field (fT^ 
{C, = 1/6, = 0), for which the de Sitter invariant Green function during inflation in D = 4 
reads j^f 

iG coni (y) = -r^- , (16) 
where y denotes the de Sitter length function 



y=-^=4sm''(^Hi], (17) 
V1V2 v2 / 

which is related to the geodesic distance i as indicated, and Ax^ = —{rji — 772)^ + ll^^^i — ^2V'- 
For an observer moving along a geodesic, y = y{x{t + At/2);x{t — At/2)) = — 4sinh^ (if At/2). 
The response function ^ for the conformal vacuum ()16|) can then be written as 

dt ^71^ J ^00 4[sinh(V2)-2e]' ' ^ ^ 

where the pole prescription corresponds to that of the Wightman function iG^. This integral 
can be easily performed by contour integration. The (double) poles (which also correspond to 
the zeros of y) all lie on the imaginary axis, u„ = ift„ = 2^7™ (ne Z). For E > Eq the contour of 
integration ought to be closed by a large circle below the real axis, such that the integral in (fTH|l 
can be evaluated by summing the residua which lie (strictly) below the real axis, as illustrated 
in figure HJ The result is 



^■^conf(AE) ^(o--. -^^E j27rAE/H)n 

dt ~ 47r2^ ^ H ^ 



_ 1 

Q{2n/H)AE _ I ' (^^) 

which is identical to the response function ^ of the thermal Bose-Einstein distribution for 
massless scalars, confirming thus the well known result [2]. For AE < 0, the contour should be 
closed above the real axis, such that the contributing poles are n > 0, also shown in figure 
The result of integration is again given by Eq. (|Tn|l . 
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FIG. 1: The integration contour for the Unruh's detector response function in conformal vacuum. The 
sohd (blue) contour corresponds to AE' > 0; the dashed (red) contour to IS.E < 0. 

B. Nearly minimally coupled light scalar 



The Green function for a massive scalar fielc 
^rangean (jl2p is given by the Chernikov-Tagirov 



coupled to gravity as indicated by the La- 
^ (Bunch-Davies ,8j) vacuum 



^G{y) 



4^ 



where 



3x2 ml + l2iH^ 
2, 



(20) 



(21) 



The uniqueness of iG{y) follows from the requirement that the lightcone singularity is of the 
Hadamard form. 

When expanded in powers of 



— V 



+ 4^ + [(mj///^) + 12^]' 



|s| < 1 , 



(22) 



2 ' 3if2 

the Green function for the Chernikov-Tagirov vacuum reduces to the following simple form 



(23) 



The nontrivial new integral comes from the term iG^^Q oc ln(?/), and its contribution to the 
response function yields the integral 
rf^in(AE) H 



dt 



In f4sinh^ + i7rsign(M) 



(24) 
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where we broke the logarithm into the real and imaginary contributions, in accordance with 
the e-prescription for iG^. The real part of the logarithm can be evaluated by breaking it into 
positive and negative u and then performing a partial integration (or, as it is more commonly 
done, by expanding the logarithm), while the imaginary part can be integrated trivially: 



dJ^U^E) r,. fAE \ ^ /n^ 

— au sm ——u coth — 



dt An'^AE Jo \ H J \2J AnAE 

1 



(25) 



27r e2'^^-^/^ - 1 ' 
where in the last step we made use of Eq. (3.981.8) of Ref. [l^ . 

The remaining integrals in the response function simply yield (5-function contributions, 

^-|?^<^^)(^-i+>"P)-«w)- 

Collecting all terms together, we get the response function for the nearly minimally coupled 
massive scalar: 



It = + A^J e(WH)A. _ 1 + ^^(^E) (- - 1 + ln(2) + 0(s^ 



(27) 



C. Minimally coupled massless scalar field 

The Green function of a minimally coupled massless scalar field exhibts an infrared diver- 
gence and the construction of a finite propagator necessarily breaks de Sitter invariance. 
For the purpose of calculating loop diagrams and using dimensional regularisation, one considers 



the propagator in d dimensions wit 



1 appropriate counterterms to cancel the infrared divergence 



For a detailed discussion, see Ref. 1^. In four dimensions, one obtains 13, 3] 

r 1 1 1 1 1 

tG^=o{xi; ^2) = ^ I - - 2 Hy) + 2 In {<ViHV2)) - ^ + ln(2) j , (28) 

where the term oc In (0(771)0(772)) breaks de Sitter invariance. 

Yet, this does not imply, that there is no de Sitter-invariant vacuum, as pointed out in 
Ref. |l2|, where such an invariant state is explicitly constructed by quantizing the mode with 
zero momentum seperately. However, singling out that mode does not render the propagator 
finite, as one sees when regulating the propagator with an infrared cutoff such that one 
obtains [3| 

^6-^=0,^0(0:1; ^2) = - ^ Hy) + ^ In {aivMv2)) - Hk'H) - 7^ + O(A;0)| , (29) 
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where je = 0.577215.. is Euler's constant. This expression differs from the propagator (j28p . 
which we shall use for calculating the response, only by a constant. 

The response function is easily reconstructed from the results of section HV Bl 

where ln(a) = Ht = N is the number of e-folds elapsed since the beginning of inflation, if we 
set the initial scale factor to be one. This contribution to the response function vanishes for all 
AE 7^ 0, such that under the assumption that its energy levels are not degenerate, the detector 
is insensitive to the breaking of de Sitter invariance by the propagator. Note also that, when 
integrated over AE around AE = 0, the terms oc 6{AE) are subdominant, because they only 
give a finite contribution to the response, provided that the scale factor a and the cutoff 
in are finite and nonzero, while the remaining terms yield a divergence. 

However, for a detector sensitive to time separations up to At < A^ax, in becomes 
a smeared "5-function," dJ^s{AE; A^s.x) / dt oc sin (AEAmax/2)/(AE/2), such that the detector 
responds roughly up to the energies AE < 27r/Amax- Let us now study the question of finite-time 
measurements in more detail. 



D. Boundary terms through finite-time measurements 

Strictly speaking, the propagators (fTI)| EBl I2HI) describe the dynamics of the scalar field for 
t > to- Hence it is natural to consider the Unruh detector, which corresponds to the transition 
amplitude Q, which begins measuring at to = and ends aXtf =t. In this case Eqs. © and (jH) 
generalize to 

— _ ii \mEQ^E\ T, {I\E - E - Eq) , 

dt dt 

with 

rfJ-(A£;,0,t) _ r^^^e-.Ai5At^^<^^^^^/2,f;t-At/2,x), 



dt 



implying that the response function gets modified by the boundary (initial) effects. For example, 
the response function associated with the Hadamard (conformal) form, iGconf = -f^^/(4vr^?/), for 
an Unruh detector reads 

dT^oAAE, 0,t) _AE 1 H ^ ncos (AEt) - {AE/H) sin (AEt) 
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Similar (though more technical) analysis can be performed for other contributions from the 
massless scalar propagator (j^HJ. Quite generically, boundary effects give rise to oscillatory con- 
tributions to the response function of Unruh's detector. For AE ~ H, these terms become 
unimportant when t ^ H^^. In the ultraviolet, where AE ^ H, the oscillatory contributions 
become subdominatnt when t ^ AE/H"^ - much more than a Hubble time. Since these effects 
are however not related to the intrinsic nature of the quantum fields in de Sitter space, we here 
do not study them further. 



E. Dimensions other than four 



So far, we have calculated the response functions from the scalar propagator, which is moti- 
vated by the practice of defining vacua in de Sitter space through this quantity. However, there is 
a method due to Higuchi [gI using a basis of wave functions as starting point, which we generalize 
here to d dimensions. We define 



d- 1 



(31) 



where the curvature is given by i? = d{d — 1)H^. The scalar wave equation is (c/. Eq. (fT^ ). 



- (1/4) 



and has the properly normalized positive frequency solution 

1 



<^k(^) 



-7rr7)^et^M/7(2)(_|k|^) 



The transition probability for the detector in terms of these modes is then 



P{AE) 



d 



(2 



TT 



OO POO 

dt / dt' 



OO J ~oo 

OO POO 



a(t)a{t 



-^tit)Mt')e- 



iAE{t'-t) 



(32) 



(33) 



(34) 



d'^-^k 
{271^- 

From this expression, one obtains for the response function 



rf-^d(AE) 



T 



{d - l)/2 + iAE/H + v\^[{d-l)/2 + iAE/H - v 



dt 87r('^+i)/2r((rf- l)/2) 

which reduces for = 4 to Higuchi's result 

'3/2 + iAE/i7 + z/ 



r 



^ H-.AE/H 

dt 47r3 



r 



r 



3/2 + iAE/H - u 



(35) 



(36) 
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For V = \j1 this coincides with our result for the conformal case (jl8|) . and when expanded in 
s = 3/2 — z/ with the nearly minimally coupled case ()27|) . 

It is however also interesting to derive some special responses from the scalar propagator, 
which is in d dimensions 

V{^^y)V{^-v) ^ (d-l , d-l d y\ 



(47r)^/2r (f) ' ' V 2 ' 2 '2' 4, 

In particular, for d = 3, the response function is exactly calculable for arbitrary u from the 
propagator, because we can express the hypergeometric function in terms of the geodesic distance 
e, using Eq. (9.121.30) of Ref. Q, as 

p ^^M ^ ^ ^ y\ sin [u (tt - H£)] 

With £ iAt, the response function (0]) is then given by the integral 

dT,{AE) _ r(z/)r(l - u)_^ I dAte-'^^^' smiTTu-tuHAt) ^^^^ 



dt 4z7r2 J smh{HAt) 

—00 

with the poles of the integrand at At„ = inn/H, ne Z\{ — 1}. Note, that for odd dimensions, the 

analytic structure of the propagator is different than for even dimensions. There are additional 

poles, but there is no branch cut (except for special values of u). We perform the integration 

by closing the contour in the lower complex half plane. According to the e-prescription for the 

Wightman function iG"^, only the poles n < — 1 contribute, and we obtain 

dJ'sjAE) ^ 1 sinh (vr^) 1__ 

dt 2 cos(7rz/) + cosh (vr^) e^-^^/^ - 1 ' ^ ' 

Note that ()4()j) applies not only when < H, and u is real, but also for > H, when 

u is imaginary. When ^ H,AE, the response is exponentially suppressed as dJ-'/dt oc 

exp(— vrm^/iir), which is a consequence of an exponential suppression of scalar particle production 

in de Sitter space in the limit when ^ H. 

Note that for no value of the parameter u the response (j4(jp agrees with the thermal response 

function, which in three dimensions can be read off from Eq. (jHj). 

% = isig„(A£)e((Ai5)^ - ,„l) . (41) 

In particular, for a conformal massless scalar, u = 1/2, and Eq. ()4()j) yields the following 



'fermionic-like' response function, 

rf^3,conf(AE) 1 1 



dt 2 e2-^^/^^ + 1 



(42) 
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This disagreement with the thermal case is not a special feature of odd dimensions. E.g. the 
conformal Green function in d dimensions 

= y^li'-V-'" («) 

leads for (i = 6 to the response 

ci^6,conf(A^) fAE\AE\ 1 



dt 127r2 \ H J e2-^^/^ + r 

while the flat-space thermal response is 



(44) 



^ = ^sign(Ai^;)0 ((AE)2 - ml) {AE^ - m'f ^ . (45) 



Generally, for d > A, the conformal response consists of the Planck factor times a polynomial 
involving different powers of AE, therefore deviating from the thermal response, which involves 
only a single power of AE. 

For d = 2 conformal and minimally massless coupled case coincide and we find from (j35|) 

C^^l+l,conf(A^) 1 1 



dt AE ef^^^ - 1 

in agreement with the flat-space thermal response. 



{AE ^ 0) , (46) 



Hence, we have found that an agreement with the thermal response occurs 
coupling only in c? = 2 and d = 4, just as for an accelerated observer in fiat space 



or conformal 



3- 



V. DETAILED BALANCE, RESPONSE FUNCTIONS AND SPECTRA 

Based on the assumption that the principle of detailed balance holds, which states that the 
absorbtion rate of the detector TZa and the emission rate Tie are equal, 

na{Eo^E)=ne{E^Eo), (VEo,E), (47) 

and on the fact that the transition probabilities per unit proper time are related as follows: 

dPjEo ^ E) ^ ^_,^^_E^^ dPiE ^ Eq) 

dt dt ' ^ ' 

or, equivalently, the response function of the detector fulfills, 

dr(^^^_,,,m_AEl f^^E^E-E„), (49) 
dt dt 
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one can infer that the detector is thermally populated, with the temperature given by T = as 
follows (for a related discussion see Ref. |4||). Let us rewrite the principle of detailed balance (jlTjl 
as 

+ = + . (50) 

where n{E) and n{Eo) denote the occupation numbers of detector states with energies E and 
Eq, respectively. From this, it immediately follows 

^(E) = _ ^ , (51) 

such that the states of the detector are populated according to a chemical equilibrium at tem- 
perature T = 1/ (3 and chemical potential /x. 
In fact, any response which can be written as 

dt " ^ 27r e^A^ - 1 ' 
with g = g{AE) being an even function of AE, fulfills the relation We have shown explicitly 
for different scalar fields in c? = 4 (c/. Eqs. p9|) . ()30p and ()27p). that they are of the form 



fi'conf — 1 



9n.,=o = 1 + j + 27rH5iAE) [Ht~- + ln(2) 



1 



1 



gm,^o = 1 + j + 27rHS{AE) - 1 + ln(2)J + 0{s) , (53) 

with Ph = ^/Th = 27!-/ H. Moreover, the more general expressions (jHEl) and pUj) also 

satisfy equation 

The relation (j49p can also be viewed as a consequence of the periodicity of the Green function 
iG"^ in imaginary proper time, t t + 27ii3 J5jl]; which is in turn a consequence of the same 
periodicity of the metric in Euclidean time j2^. An example where the periodicity of the metric 
however does not coincide with the Hawking temperature is a quasi-de Sitter space considered 
in Ref. j^], and hence cannot in general be used as an argument for the thermality of a scalar 
field. 



VI. DISCUSSION 

We found the response functions for different scalar fields to differ strongly in the infrared, 
where AE < H. Moreover, they do not in general coincide with the response to an equilibrum 
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state in flat space. A disagreement with the thermal response does not yet imply that the 
detector does not equilibrate with the de Sitter background. In fact, the energy levels of the 
detector are thermally populated. Similar deviations from a Minkowski-space thermal response 



are also known for accelerated detectors [19|, |21[. The disagreement of the response functions 
should be attributed to the fact, that for the conformally and the minimally coupled scalar field 
the density of modes per frequency is different. 

However, for fields which are massive or nonconformally coupled to the metric, the infrared 
enhancement is a consequence of the amplification of superhorizon modes leading to cosmological 
density perturbations, an effect which is absent in the conformally coupled massless case. This 
makes the different fields clearly distinguishable by observables. The fact that physically distinct 
situations such as thermally populated fiat space and the different types of de Sitter-invariant 
vacua result in the same thermal distribution function for the energy levels of the detector is due 
to the insensitivity of this quantity to whether there is a mixed (thermal) state or a pure (de 
Sitter-invariant) state and to how often the detector responds per unit time ^|. 

Finally, we want to point out that de Sitter invariance of the stress-energy tensor of the 
conformal vacuum implies that the stress-energy tensor must be proportional to the metric, such 
that Pconf = — Pconf ■ This is inconsistent with a thermal equation of state and leads us to the 
investigation, whether the energy density is captured at higher orders in perturbation theory [2^, 
allowing for the pair creation operators of the scalar field Hamiltonian, which are not captured 
by the first order perturbation expansion used here. While we argued in this paper that the 
thermal state of the detector does not contain the full information about the quantum field, it 
is yet remarkable, that at first order in perturbation theory the detector is insensitive to the 
(regulated) stress-energy tensor, which is clearly nonthermal. 



Acknowledgements 

We would like to thank Michael G. Schmidt and Richard P. Woodard for critical reading of 
the manuscript and for many suggestions, and Robert H. Brandenberger, Atsushi Higuchi, Glaus 
Kiefer, Klaus Kirsten and Bharat Ratra for useful correspondence. T.P. thanks Kiril Krasnov 
for discussions. T.P. acknowledges a generous financial support of the Albert Einstein Institute 



15 



in Golm. 



[1] W. G. Unruh, "Notes On Black Hole Evaporation," Phys. Rev. D 14 (1976) 870. 
[2] N. D. Birrell and P. C. W. Davies, "Quantum Fields In Curved Space," Cambridge University 
Press (1984). 

[3] G. W. Gibbons and S. W. Hawking, "Cosmological Event Horizons, Thermodynamics, And Particle 

Creation," Phys. Rev. D 15 (1977) 2738. 
[4] M. Spradlin, A. Strominger and A. Volovich, "Les Houches lectures on de Sitter space," published 

in "Les Houches 2001, Gravity, gauge theories and strings" 423-453 [arXiv:hep-th/0110007 ]. 
[5] R. Bousso, A. Maloney and A. Strominger, "Conformal vacua and entropy in de Sitter space," 

Phys. Rev. D 65 (2002) 104039 |arXiv:hep-th/0112218| . 
[6] A. Higuchi, "Quantization Of Scalar And Vector Fields Inside The Cosmological Event Horizon 

And Its Apphcation To Hawking Effect," Class. Quant. Grav. 4 (1987) 721. 
[7] N. A. Chernikov and E. A. Tagirov, "Quantum Theory Of Scalar Fields In De Sitter Space-Time," 

Annales Poincare Phys. Theor. A 9 (1968) 109. 
[8] T. S. Bunch and P. C. Davies, "Quantum Field Theory In De Sitter Space: Renormalization By 

Point Splitting," Proc. Roy. Soc. Lond. A 360 (1978) 117. 
[9] E. Mottola, "Particle Creation In De Sitter Space," Phys. Rev. D 31 (1985) 754. 
[10] Bruce Allen, "Vacuum States In De Sitter Space," Phys. Rev. D 32 (1985) 3136. 
[11] B. Allen and A. Folacci, "The Massless Minimally Coupled Scalar Field In De Sitter Space," Phys. 

Rev. D 35 (1987) 3771. 
[12] M. Le Bellac, "Thermal Field Theory," Cambridge University Press (1996). 

[13] T. Prokopec and E. Puchwein, "Photon mass generation during inflation: de Sitter invariant case," 

JCAP 0404 (2004) 007 a rXiv:astro-ph/0312274_ . 
[14] Izrail Solomonovich Gradshteyn, losif Moiseevich Ryzhik, Table of integrals, series, and products, 

4th edition. Academic Press, New York (1965). 
[15] V. K. Onemli and R. P. Woodard, "Quantum effects can render w j -1 on cosmological scales," 

|arXiv:gr-qc /0406098 , 

[16] T. Prokopec and R. P. Woodard, "Dynamics of super-horizon photons during inflation with vacuum 
polarization," Annals Phys. 312 (2004) 1 |arXiv:gr-qc/0310056| . 



16 



[17] K. Kirsten and J. Garriga, "Massless Minimally Coupled Fields In De Sitter Space: 0(4) Symmetric 

States Versus De Sitter Invariant Vacuum," Phys. Rev. D 48 (1993) 567 arXiv:gr-qc/9305Q13| . 
[18] N. C. Tsamis and R. P. Woodard, Class. Quant. Grav. 11 (1994) 2969. 

[19] C. Gabriel, P. Spindel, S. Massar and R. Parentani, "Interacting charged particles in an electric 



field and the Unruh effect," Phys. Rev. D 57 (1998) 6496 |arXiv:hep-th/9706030^ . 
[20] R. H. Brandenberger and R. Kahn, "Hawking Radiation In An Inflationary Universe," Phys. Lett. 
B 119 (1982) 75. 

[21] R. Brout, S. Massar, R. Parentani and P. Spindel, "A Primer for black hole quantum physics," 

Phys. Rept. 260 (1995) 329. 
[22] R. H. Brandenberger, "An Alternate Derivation Of Radiation In An Inflationary Universe," Phys. 

Lett. B 129 (1983) 397. 

[23] B. Garbrecht and T. Prokopec, "Energy density in expanding universes as seen by Unruh's detec- 



tor," |arXiv:gr- qc/0406114t 
[24] For more general^ de Sitter invariant states, the so called a- vacua, one can show that (|49)) does not 
in general hold 5]. However, these states have a different ultraviolet structure than the Chernikov- 
Tagirov vacuum H20() . which has the standard Hadamard lightcone singularity, and hence they are 
most likely unphysical. 

[25] An example for the interpretation of the de Sitter invariant states as thermal can be found in 
Ref. where scalar field quantization is performed in static coordinates. The mode functions are 
chosen to vanish beyond the horizon, where the static coordinates exhibit a coordinate singularity. 
Since the horizon distance is singled out, the mode functions in static coordinates violate spatial 
homogeneity. An Unruh detector placed at the origin of the static vacuum measures no particles. On 
the other hand, when the static vacuum is thermally populated, the response function is the same 
as for the de Sitter- invariant vacuum, and it is given by (|36() . Note first that a thermally populated 
state in static coordinate does not correspond to a usual thermal equilibrium state, since spatial 
homogeneity is broken. In addition, the mode functions in the static and the de Sitter-invariant 
vacuum have different support. Indeed, the static mode functions vanish beyond the de Sitter 
horizon, while the de Sitter-invariant mode functions exhibit superhorizon correlations. 



